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Abstract 

It has been shown recently that the normal branch of a DGP braneworld scenario self- 
accelerates if the induced gravity on the brane is modified in the spirit of f(R) modified 
gravity. Within this viewpoint, we investigate cosmological viability of the Hu-Sawicki 
type modified induced gravity. Firstly, we present a dynamical system analysis of a gen- 
eral /(.R)-DGP model. We show that in the phase space of the model, there exist three 
standard critical points; one of which is a de Sitter point corresponding to accelerating 
phase of the universe expansion. The stability of this point depends on the effective 
equation of state parameter of the curvature fluid. If we consider the curvature fluid to 
be a canonical scalar field in the equivalent scalar-tensor theory, the mentioned de Sitter 
phase is unstable, otherwise it is an attractor, stable phase. We show that the effective 
equation of state parameter of the model realizes an effective phantom-like behavior. A 
cosmographic analysis shows that this model, which admits a stable de Sitter phase in its 
expansion history, is a cosmologically viable scenario. 
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1 Introduction 



The late-time accelerating phase of the universe expansion which is supported by data related 
to the luminosity measurements of high red shift supernovae [1], measurements of degree-scale 
anisotropies in the cosmic microwave background (CMB) [2] and large scale structure (LSS) 
[3], is one of the challenging problems in the modern cosmology. The rigorous treatment of 
this phenomenon can be provided essentially in the framework of general relativity. In the 
expression of general relativity, late time acceleration can be explained either by an exotic fluid 
with large negative pressure that is dubbed as dark energy in literature, or by modifying the 
gravity itself which is dubbed as dark geometry or dark gravity proposal. The first and simplest 
candidate of dark energy is the cosmological constant, A [4]. But, there are theoretical problems 
associated with it, such as its unusual small numerical value (the fine tuning problem), no 
dynamical behavior and even its unknown origin [5]. These problems have forced cosmologists 
to introduce alternatives in which dark energy evolves during the universe evolution. Scalar field 
models with their specific features provide an interesting alternative for cosmological constant 
and can reduce the fine tuning and coincidence problems. In this respect, several candidate 
models have been proposed: quintessence scalar fields [6], phantom fields [7] and Chaplygin 
gas [8] are among these candidates. Nevertheless, we emphasize that the scalar field models of 
dark energy are not free of shortcomings. 

As an alternative for dark energy, modification of gravity can be accounted for the late 
time acceleration. Among the most popular modified gravity scenarios which may successfully 
describe the cosmic speed-up, is f{R) gravity [9,10]. Modified gravity also can be achieved 
by extra-dimensional theories in which the observable universe is a 4-dimensional brane em- 
bedded in a five-dimensional bulk. The Dvali-Gabadadze-Porrati (DGP) model is one of the 
extra-dimensional models that can describe late-time acceleration of the universe in its self- 
accelerating branch due to leakage of gravity to the extra dimension [11,12]. 

Recent observations constrain the equation of state parameter of the dark energy to be 
wx ~ —1 and even wx < —1 [13]. One of the candidates for dark energy of this kind is the 
phantom scalar field. This component has the capability to create the mentioned acceleration 
and its behavior is extremely fitted to observations. But it suffers from theoretical problems; 
it violates the null energy condition and its energy density increases with expansion of the 
universe leading to a future big rip singularity. Also it causes the quantum vacuum instabilities. 
So, some authors have attempted to realize a kind of phantom-like behavior (w e ff < — 1) 
in the cosmological models without introduction of phantom fields. In fact, the possibility 
of realization of an effective phantom nature without introduction of phantom fields is an 
important task and has been appreciated sufficiently in recent years [14]. 

In the streamline of the mentioned issues, we are going to study cosmological viability of a 
class of DGP-inspired braneworld models in which the induced gravity on the normal branch is 
modified in the spirit of f(R) gravity [10,15,16,17]. Firstly, we study the cosmological dynamics 
of this model within a dynamical system approach. We show that there exists a standard de 
Sitter point in the phase plane of the model. In this respect, this model has the potential to 
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explain accelerated expansion of the universe. The stability of this point depends completely 
on the effective equation of state parameter of the curvature fluid. If we consider the curvature 
fluid to be a canonical scalar field in the equivalent scalar-tensor theory, the mentioned de Sitter 
phase is unstable, otherwise it is an attractor, stable phase. Since the late-time accelerating 
phase of the universe expansion is explained by a stable de Sitter phase, we can investigate the 
cosmological viability of such theoretical models based on the phantom-like behavior of this 
/(.R)-DGP gravity. To be more specific, in which follows we focus on the cosmological viability 
of the Hu-Sawicki type modified induced gravity and show that this model has capability to 
realize a stable, attractor de Sitter phase. We point out that the phantom mimicry discussed in 
this study has a geometric origin. To be more realistic, we compare our results with observation 
via a cosmographic approach. 



2 /(fl)-DGP scenario 

In this section, possible modification of the induced gravity on the brane is investigated in 
the spirit of f(R) theories [10,15,16,17]. It has been shown that 4D f(R) theories in the 
present time can follow closely the expansion history of the ACDM universe [18]. Here we 
study an extension of f(R) theories to a DGP braneworld setup. The motivation behind this 
study is that modified induced gravity on the normal branch of a DGP scenario provides some 
new interesting features, one of which is self-acceleration of the normal DGP branch in this 
situation (see Refs. [10,16,17] for details). Similar to the normal branch of the standard DGP 
cosmology, the resulting generalized normal branch is also ghost-free and therefore the issue 
of ghost-instabilities is irrelevant in this case [17]. The action of this model can be written as 
follows 

ft/T 3 r. r- 

S = ^ J d 5 x^TZ + J d A x^(Ml~K + C) , (1) 

where by definition 

2 

Tfl 

C=^f(R) + C m . (2) 

By calculating the bulk-brane Einstein's equations and using a spatially flat FRW line element, 
the following modified Friedmann equation is obtained [15,16,17] 



H2= 8nG / (m) + p{rad) + ( a ± H 
3 V / r c 



= 87rGY (m) | Jrad) | j curv) \ , H 

where 



p (curv) = m 2 



f(R) - Rf'(R)} - 3RHf"(R)j , (4) 



1 

is energy density corresponding to the curvature part of the theory. This energy density can be 
dubbed as dark curvature energy density. f c is the re-scaled crossover distance that is defined 
as f c = r c f'(R) and a prime marks differentiation with respect to the Ricci scalar, R. We note 
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that in this scenario there is an effective gravitational constant, which is re-scaled by f'(R) so 
that G = G e ff = 8nm 2fi( R ) [15]. In order to study the phase space of this scenario, it is more 
suitable to rewrite the normal branch of the Friedmann equation (3) in the following more 
phenomenological form 



e 2 = n m (i + z f + n rad (i + z f + n curv (i + z f^+^ v ) _ 2 ^n rc E , (5) 

where by definition 

_ 8ttG (curv) I 

J L curv j_ T o HO ) J L r 



and also 



Rf"(R) + R Rf"'(R) - Hf'(R) 



Uf(R)-Rf'(R)]-3HRf»(R) 



(6) 
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We note that w curv is not a constant and varies with redshift. 



3 The phase space of a general /(i?)-DGP model 

To investigate cosmological dynamics of this model within a dynamical system approach, we 
express the cosmological equations in the form of an autonomous, dynamical system. For this 
purpose, we define the following normalized expansion variables 



V 



a 3 / 2 E ' 



r = 



rad 



a 2 E 



s = 



a 3{l+w curv )/2 ft ' 



U = 



E 



(7) 



In this way, equation (5) with minus sign (corresponding to the generalized normal DGP branch) 
and in a dimensionless form, is written as follows 



1 + 2u = p 2 + r 2 + s 2 . 



(8) 



This constraint means that the allowable phase space of this scenario in the p-r-s space is 
outside of a sphere with radius 1, which is defined as p 2 + r 2 + s 2 > 1 . The autonomous system 
is obtained as follows 



p 



3p 



p 2 +(l+2w cur v)r 2 



2( p 2 +r 2 + s 2 + 1 ) 



3r 2p 2 + h 2 + (l+w curv )(r 2 -p 2 -s 2 -l) 



s' = 



2{p 2 +r 2 +s 2 +l) 



2s 2 +p 2 +(l+3w CU rv)r 2 -2 



(9) 



p 2 +r 2 +s 2 +l 
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Table 1: Eigenvalues and the stability properties of the critical points. 



points (p, r, s) character 



eigenvalues 



^ curv ^ 1 ^ curv ^ 1 



A (0, 0, 1) radiation 
B (1,0,0) matter 
C (0, 1, 0) de Sitter 



1 _3 3 

2 1 2 curv 1 2 



) 3(l+4w eIlr „) 
2 ' 2 ' 8 



unstable unstable 
unstable unstable 
stable unstable 



Here a prime marks differentiation with respect to the new time variable r = In a that a is 
scale factor of the universe. The critical points in the phase plane are obtained by solving the 
equations p' — 0, r' — and s' = 0, that is, setting the autonomous system (9) to be vanishing. 
The results are shown in table 1. To investigate the stability of these points, we apply the linear 
approximation analysis to achieve the Jacobian matrix. Note that the critical points and their 
stability depend on the value of w curv . Here we investigate the stability of the standard critical 
points in two different subspaces of the model parameter space where EoS of the curvature fluid 
has either a phantom or a quintessence character. As we see in table 1, the radiation dominated 
phase (point A) and matter dominated phase (point B) in this scenario, are unstable phases of 
the universe expansion independent on the value of w curv . Whereas, the accelerating phase of 
the universe expansion (point C) is a stable phase if the curvature fluid is considered to be a 
non-canonical (phantom) scalar field (w curv < —1) in equivalent scalar-tensor theory; otherwise 
it is an unstable phase. It is necessary to mention that whenever w curv = —1, the variables 
s and u are not independent and the phase space is 2D (here the curvature fluid plays the 
role of a cosmological constant, the same as ADGP model. For more details see Ref. [19]). 
Figure 1 shows the 3D phase space trajectories of the model. In this figure, the point C as a 
de Sitter point is an attractor for w curv < — 1. Therefore, a model universe which is described 
by modified induced gravity on the normal DGP branch, has a stable, positively accelerated 
expansion phase if the modified gravity indicates a phantom-like behavior. We note that points 
O and C do not belong to physical phase space of our model universe. 

After exploration of the cosmological dynamics in a general /(-R)-DGP setup within a phase 
space analysis, in the next section we study cosmological viability of an specific f(R)-DGP 
model. 
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Figure 1: The 3D phase space of the autonomous system (9) for w curv < — 1. There are three critical 
points: C is an attractor de Sitter point, O is a saddle point and C is an attractor point. 

4 Cosmological viability of the Hu-Sawicki type modi- 
fied induced gravity 

Now we focuss on the cosmological viability of the model by considering a Hu-Sawicki type 
modified induced gravity on the DGP brane. It is shown in the Ref. [18] that the expansion 
history of the mentioned model in 4 dimensions is widely close to the ACDM model in the 
high-redshift regime. Now in a braneworld extension, we expect the Hu-Sawicki induced gravity 
mimics the ADGP model in the mentioned regime. In other words, in this regime curvature 
term in the Friedmann equation is close to the cosmological constant which is screened by the 
term ¥- . In fact, the dynamical screening effect is the main origin of the phantom-like behavior 
of the curvature term in the normal branch of this DGP-inspired braneworld scenario [15]. The 
Hu-Sawicki model [18] is given by 

f(R) = R-m 2 , (10) 

where m 2 , c\, C2 and n are free positive parameters that can be expressed as functions of density 
parameters. Now we explore the dependence of these parameters on density parameters defined 
in our setup. Variation of the action (1) with respect to the metric yields the induced modified 
Einstein equations on the brane 

G a/ 3 = J^SaP - , (H) 

where £ a p (which we neglect it in our forthcoming arguments), is the projection of the bulk 
Weyl tensor on the brane 
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and S a/ 3 as the quadratic energy-momentum correction into Einstein field equations is defined 
as follows 

111 1 

<S a /3 = V| + —TT a/3 + -gapT^T^ - — g af3 T 2 . (13) 

r a/ 3 as the effective energy-momentum tensor localized on the brane is defined as [10] 



771 r i r n 

r af s = -m 2 J\R)G afi + -f[f(R) - Rf'(R)]g a p + T a p + m 2 p [V a V p f{R) - g aP Uf\R)\ . (14) 



The trace of Eq. (11), which can be interpreted as the equation of motion for f'(R) , is obtained 

as 



R = 



24M| 



(15) 



r, the trace of the effective energy-momentum tensor localized on the brane is expressed as 



T = 771, 



2f(R) - Rf'(R) - 3Df\R) - (p m + p rad ) , 



(16) 



To highlight the DGP character of this generalized setup, we express the results in terms of the 

2 

DGP crossover scale defined as r c = ^fs- So, the equation of motion for f'(R) is rewritten as 
follows 

\r 2 ([2f{R) - Rf'(R)] 2 + 9(o f(R)) 2 + QRf'(R)Df'(R) - 12f(R)Df'(R)\ 



+ 



6 v 

5 r c 
6M, 3 



Rf'(R) - 2f(R) + 3nf(R) (p m + p ra 



24M| 



{Pm + Pradf ~ R = 0. 



(17) 



In the next stage, we solve this equation for Df'(R) to obtain 



*f{R) = -3 



(Rf'(R) - 2f(R)) + 



Prn ~T Prad 

2Mfr c 



^vV "15 ' 



(18) 



Now we introduce an effective potential V e ff which satisfies Of'(R) = §jr^ ■ This effective 
potential has an extremum at 



Rf'(R) - 2f(R) 



1 1/6 

— {Pm + Prad) = ± — \\-R- 



mi 




r, V 5 



(19) 



In the high-curvature regime, where f'(R) — 1 and ~ 1 , we recover the standard DGP 
result (one can compare this result with corresponding result obtained in Ref. [18] to see the 
differences in this extended braneworld scenario) 



1/6 1 

R±-^ l R=^{ Pm + Prad). 



(20) 



The negative and positive signs in this equation are corresponding to the DGP self-accelerating 
and normal branches respectively. In which follows, we adopt the positive sign corresponding 
to the normal branch of the scenario. To investigate the expansion history of the universe in 
this setup, we restrict ourselves to those values of the model parameters that yield expansion 
histories which are observationally viable. We note that the Hu-Sawicki f(R) function intro- 
duced in Ref. [18], was interpreted as a cosmological constant in the high-curvature regime. 
The motivation for that interpretation was to obtain a ACDM behavior in the high curvature 
(in comparison with m 2 ) regime. Here we show that in a braneworld extension, the Hu-Sawicki 
induced gravity mimics the ADGP model in the mentioned regime. As we have pointed out 
previously, the phantom-like behavior can be realized from the dynamical screening of the brane 
cosmological constant. In this respect, we apply the same strategy to our model, so that the 
second term in the Hu-Sawicki f(R) function (that is, the second term in the right hand side of 
equation (10)) mimics the role of an effective cosmological constant on the DGP brane. Then 
this term will be screened by =- term in the late time (see the normal branch of Eq. (3)). 
In the case in which R ^> m? , one can approximate Eq. (10) as follows 

lim f(R)^R-^ m 2 + % m 2 (^-) n . (21) 

z£_> c 2 c| 

During the late-time acceleration epoch, f^R) — 1 or equivalently R S> m 2 and we can apply 
the above approximation. Also the curvature field is always near the minimum of the effective 
potential. So, based on Eq. (19), we have 

= \ {pm + Prad) + 2-m 2 . (22) 

Since R in the f(R) function is induced Ricci scalar on the brane, we except crossover scale 
to affect on the constant parameters C\ , c 2 and m 2 . In Ref. [18] the authors obtained 3m 2 = 
Rc = that pom is the present value of the matter density. But, in our setup the present 

value of the matter density (see Eq. (20)) is given by 

[R~ 1 

R c + 0.9^ = — ( POm + POrad) ■ (23) 

r c m 2 p 

If we solve this equation for R c , we find 

3m 2 = R c « Q rc + 3Q m + 3Q rad ± 3^ rc (o.068fi rc + Q m + Q rad ) . (24) 

Therefore, the DGP character of this extended modified gravity scenario is addressed through 
m 2 . As we have mentioned, at the curvatures high compared with m 2 , the second term on the 
right hand side of equation (10) mimics the role of an effective cosmological constant on the 
brane. In this respect, the second term in the right hand side of equation (21) also mimics the 
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role of a cosmological constant on the brane in the high curvature regime. With this motivation, 
we find 

|« 18 " A (25) 



There is also a relation for % as follows 

ci _ l-^(i?) /i? y+ 1 



(26) 



where ^ in our setup can be calculated as follows: firstly, by using Eqs. (22) and (25), we find 



1 6 „ 1 



R+ -\h R= — (Pom a + POrad O + 12fi A , (27) 



r r V 5 m 2 



where p om can be omitted through Eq. (23) to obtain 



R + -\I^R = ( 3m 2 + 1.56- - — p 0rad )a~ 3 + — p 0rad a~ 4 + 12fi A . (28) 
rJ 5 V r r mi J mi 



l/6„ /„o . .„rn 1 \ _ 3 1 

~Z POrad O H 

Finally, if we solve this equation for V^R , we find the following relation for ^ 



R _ ( _ 0.9^ + 



m 2 \ m 



, | 0.52^ x2 | V2Q,l 

V m ' 



m ' m 2 



2 



(29) 



where m is given by Eq. (24). Note that we have set Hq and a(to) equal to unity. These 
relations tell us that the free parameters of this model are n, Q m , Q r ad, ^r c and fb(R), whereas 
the latter one is constrained by the Solar-System tests. In fact, experimental data show that 
f'(R) — 1 < 1CT 6 , when f'(R) is parameterized to be exactly 1 in the far past. To analyze 
the behavior of w curv , we specify the following ansatz for the scale factor 

a w = (' 2+ r^r- (30) 

where v ^ 1 is a free parameter [20]. By noting that the Ricci scalar is R = 6(- + (f ) 2 ), one 
can express the function f(R) of equation (10) in terms of the redshift z. Figure 2 shows the 
variation of the effective equation of state parameter versus the redshift. As we see in this 
figure, in this class of models the curvature fluid has an effective phantom-like equation of 
state, w curv < —1, at high redshifts and then approaches the phantom divide {w curv = —1) at 
a redshift that increases by decreasing n. 

The main point here is that a modified induced gravity of the Hu-Sawicki type in the DGP 
framework, gives an effective phantom-like equation of state parameter for all values of n, and 
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Figure 2: w CU rv versus the redshift for a Hu-Sawicki type modified induced gravity with O m = 0.27, 
&rad = 0.3, $7a = 0.9, fl rc = 0.01 and fb(R) — 1 ~ 10 -6 . As this figure shows, in this class of models 
the curvature fluid has an effective phantom equation of state with w curv < — 1 at high redshifts. This 
effective equation of state parameter approaches the phantom divide line (w curv = — 1) at a redshift 
that increases by decreasing the values of n. 

all of these models approach asymptotically to the de Sitter phase {w curv = —1). Therefore, 
the accelerated expansion of the universe (the de Sitter phase) is necessarily a stable attractor 
phase for this DGP-inspired f(R) model. Based on the analysis presented in the previous section 
within a phase space viewpoint and also the outcomes of this section, we can conclude that a 
Hu-Sawicki type modified induced gravity on the normal branch of the DGP setup provides 
a cosmologically viable scenario. This is the case since it contains a radiation dominated era 
followed by a matter dominated era and finally a stable de Sitter phase in its expansion history. 
In the next section we compare our model with observational data via a cosmographic analysis. 
Our treatment here is mainly based on the Ref. [29,30]. 

5 Comparison with observational data 

While theoretical consistency of a physical theory is a primary condition for viability of the 
theory, the observational consistency of the model is necessary too. For this goal, in which 
follows we discuss briefly observational status of our model via a cosmographic analysis. Before 
that, we note that the DGP model is a testable scenario with the same number of parameters 
as the standard ACDM model, and has been constrained from many observational data, such 
as the SNe la data set [21], the baryon mass fraction in clusters of galaxies from the X-ray 
gas observation [22], CMB data [23], the large scale structures [24] and the baryon acoustic 
oscillation (BAO) peak [25], the observed Hubble parameter H(z) data [26], the gravitational 
lensing surveys [27]. The observational constraints on the DGP model with Gamma-ray bursts 
(GRBs) at high redshift also obtained recently from the Union2 Type la supernovae data 
set [28]. In [28] the authors are shown that by combining the GRBs at high redshift with 
the Union2 data set, the WMAP7 results, the BAO observation, the clusters baryon mass 
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fraction, and the observed Hubble parameter data set and also in order to favor a flat universe, 
the best fit of the density parameter values of the DGP model are obtained as {fi m , Q r } = 

{o.235±8:8iS ,o.i38±8:8S} [28]. 

Here to compare our /(i?)-DGP model with observational data we adopt the cosmography 
approach. Cosmography approach is a useful tool in order to constrain higher order gravity ob- 
servationally without need to solve field equations or addressing complicated problems related 
to the growth of perturbations [29,30]. In this case, one can define cosmographic parameters 
based on the fifth order Taylor expansion of the scale factor. One can also relate the charac- 
teristic quantities defining the f(R)-DGP model to the mentioned cosmographic parameters. 
Therefore, a measurement of the cosmographic parameters makes it possible to put constraints 
on f(Ro) and its first three derivatives. The likelihood function for the probe s is defined as 
[31] 

£ s (p)ocexp(- X 2(p)/ 2 ) (31) 



where 



2 



X 2 S (P) = E 1 — (32) 

71=1 U ^\^l) 

Hobs(zi) are the observed distance modulus for the adopted standard candle (such as SNe la) at 
the redshift Zi with its error o^. fith{z) are the theoretical values of the distance modulus from 
cosmological models which read as /i t h(z,p) = 25 + 5 log-D^z, p) where Dl = H cIl is the 
luminosity distance. In the cosmography approach, one can obtain an analytical expression for 
luminosity distance versus the cosmographic parameters so that one require no priori model to 
solve d,L — (1 + z) J Q -gT^Ty . By using the least squares fitting that means the getting of xtmim 
one can obtain the suitable cosmographic parameters. In the next step, one should relate the 
f(R) function and its first three derivatives to the cosmographic parameters to set constraints 
on the parameters of the f(R) function [29,30]. In this manner we constrain observationally 
the parameters of a Hu-Sawicki type f(R) induced gravity on the normal DGP brane by the 
cosmography approach. Our strategy in this cosmographic approach is mainly based on the 
recent paper by Bouhmadi-Lopez et al. [30]. Firstly we relate the functions f(Ro), f'(Ro), 
f"(Ro) and f'"(R ) to the parameters R , R , R , (^r)o and H which are expressed versus 
the cosmographic parameters by using the Friedmann and Raychaudhuri equations at t — to . 
Now we have a system of two equations with four unknowns. To expand the f(R) function and 
its derivatives versus these cosmographic parameters, we need to two further equations to close 
the system. In 4-dimensional f(R) gravity, the Newtonian gravitational constant G is replaced 
by an effective (time dependent) quantity as G e // — ~WrS- On the other hand, it is reasonable 
to assume that the present day value of G e // is the same as the Newtonian one G e ff(z = 0) = G 
or f'(Ro) — 1 • One may note that the Hu-Sawicki model with this condition reduces to the 
Einstein-Hilbert gravity with Lagrangian f(R) = R. In order to resolve this problem, we can 
replace the condition f'(Ro) = 1 with f'(Ro) = (l + e) _1 . Another relation can also be obtained 
by differentiating the Raychaudhuri equation [29,30]. We solve this system of four equations 



11 



for four unknowns to obtain the following relations 



/Cgo) _ 
f"'(Ro) . 

W%r 2 ' 



V 



+ e(2 - g ) , 



A 2 n m +B 2 +C2(r c Ho)- 1 C B 2 -C 2 



V 



(j -q -2)V 



+ e 



v 



C/o-go-2)x> 



(33) 



where «4j, Z^, Cj and D with i = 0,2,3 are functions of go, jo, s o and Iq (these functions are 
defined in Ref. [30]). The new quantities C 2 , B' 3 , and C' 3 are defined as follows 



C 2 = jo(jo - <?o - 1) + Qoilo + go - 3) - 2 , 



(34) 



(35) 



B' 3 = 2j (2g 2 + 6g + j + 3) + 2g (15g 2 + 42g + 39) 
-2Zo(l + go)-2so(g +jo) + 24, 

C 3 = jo( - Jo + [2g + 8]g + s + 7) + s (l - g 2 ) 

(go[g 2 + 8g - 2] - s - 17) + 8 . (36) 

In the second step we have to determine the values of the cosmographic parameters that have 
the best fit to the observational data (by the least squares fitting). Instead, here we use 
a minimal approach to parameterize the cosmographic parameters by the phenomenological 
density parameters. In other words, the cosmographic parameters will be calculated for a given 
phenomenologically parameterized dark energy model. The best choice is the ACDM model. In 
Ref. [30] the details of these calculations are done. They finally obtained the following results 
[30] 

3 9 27 

g = -l + -^ m , Jo = l, s = l--tt m , l = 1 + 3fi m + — . (37) 



Now one can substitute these results into equations (33) and consider the observational conser- 
vative values Q m = 0.266 and f2 rc = 10~ 4 where f2 rc = (ArlHl)^ 1 [2]. Finally, by considering 
the first order approximation in e, one obtains the following results 



/(go) 



0.849 + 2.6e, 



f"{Ro 
(6# 2 )- 



0.16-20.5e, 



f"{Ro 



1.3 + 0.0176e 



(38) 



In the HS model, there are four parameters Ci, C2, R c and n that can be constrained by 
observational data via the values of the f(Ro) and its derivatives. So, we should create a system 
of four equations in the four unknowns through equation (10) and its first three derivatives. 
By solving equation (10) and its first derivative for c± and c 2 , with f(Ro) = 0.849 + 2.6e and 
f'(R ) = (1 + e)- 1 , one finds 



Ci = c x R 



l-n 
c 



n(l + e) pl _ n 
K o 



1 - 



0.849 + 2.6e 



R 



(39) 
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c 2 = c 2 R: 



n(l + e) 



D — 71 

K 



1 - 



0.849 + 2.6e 



(40) 



i?o n(l + e). 

By substituting relations (39) and (40) in HS f(R) function and its derivatives, it is obvious 
that parameter R c cancels out so that we have to work with two parameters Z\ and c 2 instead 
of three parameters ci, c 2 and R c . In other words, R c cannot be obtained in this fashion. By 
setting the second derivative of the HS function equal to f"(Ro) = 0.16 — 0.5e, we get 



[(0.849 + 2.Qe)/R ] + [(1 + e)/e](l - [0.16 - 0.5e]/i? ) - (1 - e)/(l + e) 

ft —— 

1 - (0.849 + 2.6e)/i? 



(41) 



In the last stage and in order to determine the value of e , one can use the third derivative of 
the HS function and setting f"'(R ) = 1.3 + 0.0176e to obtain the following constraint (see also 
[29]) 



1.3 + 0.0176e _ l + e(0.16-0.5e) 
0.849 + 2.6e ~ ~~ e 



Ro 



Rq 



(0.16 -0.5e) e(0.849 + 2.6e) 
0.849 + 2.6e +_ ~~ 



2e 



1 + e 



V 1 - (0.849 + 2.Qe)/R 

(42) 

Using this constrain, the acceptable value of e is ~ 0.03 (note that there are three other 
values of e that are not acceptable since are very large). The value of R is determined by 
Rq = 6Hq(1 + q ) with q = —1 + §f2 m . By equation (41), we get n ~ 2 and by equations 
(39) and (40) we obtain Z\ ~ 10 and c 2 — 0.7. Note that as we excepted these parameters are 
positive. The parameter R c here plays the role of a scaling parameter. We obtain c\ and c 2 
from equations (25) and (26) and then by using their relation with C\ and c 2 , we find R c ~ 0.018 
which is a reasonable value for this quantity. 



6 Summary 

Recently a mechanism to self-accelerate the normal branch of the DGP model, which is known 
to be free from the ghost instabilities, has been reported [17]. This mechanism is based on 
the modified induced gravity. In this paper, firstly we studied the cosmological dynamics of 
this model within a phase space approach. A de Sitter phase is the simplest cosmological 
solution that exhibits acceleration. As we have shown in a dynamical system viewpoint, this 
phase appears in our generalized setup. In fact, based on the dynamical system approach, we 
showed that there exists a de Sitter fixed point in phase space of a general /(.R)-DGP model. 
In order to investigate the stability of this accelerating phase of expansion, we classified the 
f(R) functions in two different subspaces of the model parameter space. We have shown that 
if the f(R) induced gravity plays effectively the role of a phantom scalar field in the equivalent 
scalar-tensor theory, it leads to a stable de Sitter solution and these models are cosmologically 
viable. Then, as an specific model, we studied the Hu-Sawicki type modified induced gravity 
in the DGP framework and we found that the equation of state parameter of the curvature 
fluid has an effective phantom-like character. The origin of the phantom-like behavior in the 
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model presented here can be due to the dynamical screening effect of the curvature term (which 
plays effectively the role of a cosmological constant in high-redshift regime on the brane). In 
other words, in this case the phantom-like behavior has a pure gravitational origin. We have 
shown also that the Hu-Sawicki modified induced gravity mimics the ADGP model in the 
high-redshift regime. Since the Hu-Sawicki modified induced gravity contains an early time 
radiation dominated era followed by a matter domination era and then a stable de Sitter phase 
in its expansion history, it is cosmologically a viable scenario. This result is independent on 
the value of free parameter n of the Hu-Sawicki model. Finally we have tried to constrain our 
model based on the observational data through a cosmographic procedure. In this manner we 
obtained reasonable values for parameters of the Hu-Sawicki induced gravity. 
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